Current-voltage relation for superconducting d-wave junctions 
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We calculate the current- voltage (I — V) relation for planar 
superconducting d-wave junctions for both arbitrary trans- 
mission of the junction and arbitrary orientation of the d-wave 
superconductors. The midgap states (MGS) present at inter- 
faces/surfaces of a d-wave superconductor influence the I — V 
relation. In some arrangements we find considerable nega- 
tive differential conductance and lower threshold voltage for 
non-zero current due to resonant conduction through MGS. 

PACS numbers: 74.50.+r, 74.25, 74.20.-z 



Soon after the discovery of the high-T c cuprates it was 
realized that these materials did not in a straightforward 
way conform to the conventional BCS scheme. Therefore, 
a number of researchers started to speculate whether the 
symmetry of the order parameter in the high-T c cuprates 
might deviate from the spherical symmetric one (s-wave) 
found in the conventional superconductors. Recently, an 
already lively discussion on the subject has become even 
more intensified. Quite a few experiments probing the 
magnitude of the order parameters do at least not ex- 
clude the possibility of an order parameter with symme- 
try different from s-wave. In fact, some experiments sen- 
sitive to the phase of the order parameter are claimed to 
be impossible to understand without the..assumption of 
d-wave symmetry of the order parameter raj On the other 
hand, there are papers claiming that there is a measur- 
able (if not, dominating) s-wave component of the order 
parameter aa Besides the experimental activity, there are 
several microscopic theories formulated that advocate d- 
wave symmetry, in part using arguments based on Jjhe 
strong magnetic interactions present in the cuprates.^! 

In addition to the discussion above, Hu has pointed 
out a possible connection between the often observed 
zero-bias conductance peaks in high-T c junctions and the 
zero-energy bound states or midgap states (MGS) theo- 
retically predicted to be .present at interfaces/surfaces of 
d-wave superconductorsJ3 The origin of MGS is due to 
normal scattering at the interface/surface. A quasipar- 
ticle in the superconductor changes its momentum when 
scattered. Therefore, the quasiparticle in general expe- 
riences a different order parameter before and after the 
scattering event, since in the d-wave case the order pa- 
rameter depends on the momentum. Due to Andreev 
reflection the quasiparticle will retrace its path and form 
a closed loop, leading to a bound state when the order 



parameter before and after scattering differs in sign.l 
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FIG. 1. Layout of the junction in the afo-plane. The 
d-wave superconductors "1" and "2" can be rotated an an- 
gle ai.,2. A quasiparticle is incident on the barrier at an angle 
0. The length L is only conceptual and set to zero in this 
paper. 

Calculations of I — V curves for normal metal- 
superconductor (NS) junctions indeed_cpnfirmed the im- 
portance of MGS in the d-wave case.0i3 Due to conduc- 
tion through MGS, a peak at zero voltage grows up in 
the conductance as the transmission of the junction is 
decreased, which is different from the s-wave case. 

Motivated by both the general discussion regarding 
the underlying symmetry of the order parameter of the 
cuprates and the new understanding of MGS in d-wave 
junctions we calculate in this paper the current-voltage 
(I — V) relation for junctions with anisotropic supercon- 
ductors. To produce numerical results allowing us to 
compare our calculation to earlier calculations we special- 
ize to the case of d-wave symmetry. Therefore we choose 
the order parameter to be A(6>) = A cos(2(#— a)), where 
9 is the angle of incidence for a quasiparticle approach- 
ing the junction and a is the angle of orientation of the 
d-wave superconductor with respect to the interface as 
explained in Fig. Yu. The cases a = and a = ir/4 corre- 
spond to d x 2_ y 2 and d xy symmetry, respectively. For the 
s-wave case, A(0) = A s . 

There are some calculations_of__[ — V curves for d- 
wave junctions in the litpjatureOlij Using the tunneling 
hamiltonian formalism, tM3 one can calculate the con- 
ductance as a convolution of the density of states in the 
d-wave superconductors. Unmodified, such an approach 
cannot describe the action of MGS. Now, in some con- 



figurations of <i-wave superconductors there are no MGS, 
and we can therefore compare the curves with our cal- 
culations (in the tunneling limit) as we will do later in 
this paper. Refs. [lj and [lj are using methods somewhat 
similar to ours; however, the presence of MGS is not 
discussed. There is also work on asymmetric (different 
gaps of the electrodes) SNS junctions in point-contact 
geomet™. lone or few conducting modes in the normal 
region) JiMa This theory applies straight to the d-wave 
case since for point-contacts the modes making it through 
the contact correspond to 6 » 0. In this case the MGS 
are absent. However, for the planar junction quasiparti- 
cles with any angle 9 of incidence come into play, making 
it necessary to consider the MGS. 

When there are superconductors on both sides of the 
junctkip- multiple Andreev reflection (MAR) will take 
place.cEj In this paper we use methods recently devel- 
oped to calculate I — V curves of SNS^uipctions for arbi- 
trary transmission in the s-wave caseflEHl These methods 
will automatically take into account any effects from the 
MGS, if appropriately modified. The implementation of 
these modifications is our new contribution, as presented 
in this paper. 

Our calculation is two-dimensional, using a cylindrical 
Fermi surface .0 Furthermore, the normal scattering will 
change the angle (or momentum) of an quasiparticle from 
9 to 9 = 7r — 9. This simply means that the quasiparticle 
experiences a different pair potential before and after the 
scattering event, since in general A(0) ^ A(0). (The 
exception is the case a = 0.) Therefore, one has .to book- 
keep the two different gaps in the calculation!] tl 

The system we study is shown in Fig. |l| The regions 
"Left" and "Right", separated by a barrier with trans- 
mission amplitude t(9) and reflection amplitude r(9), are 
only conceptual since we will only study the case when 
the normal region is short (zero- length limit). We model 
the barrier using a (5-functkm-potential to derive the an- 
gular dependence of r and i.113 

The time-dependent Bogoliubov-de Gennes equation 
for anisotropic superconductors is 

f fco(x,y) A(x,y,t) } ( u(y,t) ) _ 
dy { A*(x,y,t) -fto(x,y) J U(y,i) J " 

^Cm)' Mx,y)=*(x-y) ( f|-M). (1) 

In the equation above, /i is the chemical potential, which 
is [i\ {\x%) in superconductor 1 (2). Since there is a volt- 
age eV between the two superconductors, the chemical 
potentials are different , eV = fj,2 — Hi- Due to this differ- 
ence, the phase <f> of the order parameter | Ale^ depends 
on time according to the Josephson relation dt<f> — 2eV/Ti, 
leading to inelastic scattering. (For details, see Ref. Eq .) 
Solving Eq. (|l|) piecewise in each region, we follow 
the standard treatment previously used in a number of 
papersflia To this end, the center of mass coordinate 
R = (x + y)/2 and the relative coordinate r = x — y are 
introduced, expressing the order parameter as A(r, R, t). 



We will in the following assume approximate solutions 
to Eq. (jl]) of the form (in each region) 



( u(x,t) 1 
I w(x,t) J 
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iE k t 
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We approximate the integral operator of Eq. ([j]) by ne- 
glecting terms o£-.the order (feF^o) -1 (the quasiclassical 
approximation) £3 



f dyA(x, y)v(y, t)»^ A(k, x)v k e' 



k-x e -2E k t 



(3) 



where we have introduced the Fourier transform A(k, R) 
of the order parameter A(r, R). Since in the weak- 
coupling limit the pair potential is expected to be non- 
zero only close to kp, one can replace the momentum 
by an angle 9. Neglecting self-consistency, we therefore 
model the order parameter as 



Ai(0) 
A(k,x,t) = { 0, 

A 2 (9)e 



i{<j>0+2eVt/h) 



X < -L/2 

\x\ < L/2 (4) 
x > L/2. 



Since the overall phase is arbitrary, Ai is chosen real. 

Before solving Eq. (pi), it is convenient to introduce 
some notational simplifications. First we define the BCS 
coherence factors 



Vi(E, 9) E- sgn(E)y/E*-Ai(0)* 



Ui (E,6) 
Vi{E,9) 
Ul (E,9) 



A;(0) 



\E\ > W(0)\ 



E-i^A t {9) 2 -E 2 
M9) 



, \E\ < \\(9)\ (5) 



where \ui\ 2 + \vi\ 2 = 1 and i = 1,2 refers to super- 
conductor 1 and 2. Using this definition, one can ex- 
press the Andreev reflection amplitude at superconduc- 
tor i as Ai(E n ,9) = A^ n = v i {E ni 9)/u l (E ni 9) and 
Ai(E n ,9) = A i}n , where E n = E + neV. The transmis- 
sion amplitude for an electron-like quasiparticle at en- 
ergy E and angle 9 from superconductor i to enter the 
electron branch in the normal region is Ji(E,9) = Ji = 
(u 2 (E,9)-vf(E,9))/u t (E,9). 

Solving the scattering problem, we follow the ap- 
proach by Averin and Bardas.li3 First, we assume that 
an electron-like quasiparticle is approaching the barrier 
from the left superconductor (labeled by — >) at energy E 
and angle 9. Since we will calculate the current in the 
normal region to the left of the barrier, see Fig. n], we 
need the wave function ^>l(E,9) = ^>l that solves Eq. 
(0) in region "Left" (L), see Fig. @: 
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where k = (k x ,k y ) = k(cos9,siii9) 1 k = (~k x ,k y ) — 
k(cos9,sm9) and E n = E + neV (n is an even integer). 



Matching the wave functions of the various regions at 
the two NS interfaces and the barrier, the following (re- 
cursive) relations can be derived that determine the co- 
efficients a, b, c and d. For d we have 

Otndn+2 + 0nd„ + Jndn-2 = rJi8 n ,0, (7) 

where ot n , 0n an d In are defined as 

,2 ^2,n+l^-l,n+2 
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The coefficient & is determined from 

\t\ 2 A 2 , n +ld n + Ai in+2 (|r| 2 - A 2 , n+ lA 2 ,n+l)d n+ 2 



b„ 



+2 



r(l - A 2 ,„+iA 2 ,„+i) 



(9) 



The coefficients a and c are related to b and d through 

o n = Ai. n b n + Ji(5„,o, c„ = Ai )Tl d„ . (10) 

Using the wave function in Eq. (o) we can now calcu- 
late the current / per a6-plane (repeating the calculation 
above for the left-movers): 
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(11) 



where T = J d9\t\ 2 cos 0/2 and L y is the junction length 
in the ^-direction. The quantities T T (E, 9) are defined as 

T T (£, 0) = /V T (£, 0) [f(E)T;(E, 9) - f{-E)T£{E, 9)] , 
f(E)=l/[l + exp(E/k B T)], (12) 

where N T (E, 9) is the bulk density of states evaluated in 
superconductor 1 (2) for r =— ► (<— ), and 



Tj(s,e) = ^[|<j 2 -K| 2 ], 

n 

ijr(^«)=y;[i6;i 2 



„T|21 



(13) 



The Eqs. (|7|)-(|l0[) are the main technical results of this 
paper. Together with Eq. (O) we will in the following use 
these equations to numerically calculate the I—V relation 
for some d-wave junctions. Putting A = A (valid for s- 
wave) and assuming Ai = A 2 , Eqs. (|7|)-(jl0|) will exactly 
transform to Eq. (5) of Ref. [1(1 
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FIG. 2. Varying the transmission T, we show I — V curves 
for some arrangements of s-wave and d-wave superconductors. 
In (a) we show the results for the N/d xy junction. In (b)-(f), 
the I — V curves for junctions with superconducting electrodes 
on both sides are shown: in (b) s/d xy with A s = O.lAo; in (c) 
s/d xy with A s = A ; in (d) d x 2_ y 2 /d x 2_ y 2\ in (e) d xy /d xy ; in 
(f) d x 2 _ y 2 / d xy . Zero-temperature is assumed. 

In Fig. D we present some calculations of J — V curves 
involving d- wave superconductors. The MGS shows up 
in Eq. (|S|). The quantity 1 — A n A n found in the de- 
nominators of Eq. (|J) is zero when E n — for the case 
A(0) = —A(9). This leads to transmission resonances 
which we explore in the following. 

For the ballistic curves (transmission T = 1, no MGS 
in this case) involving two superconductors there is a non- 
vansishing zero- voltage curxepl that disappears when the 
transmission is not unity.EJ'E^I In Fig. 0(a) we confirm 
that using the Eos. p)-(pT[) reproduces previous results 
for the NS case.LrEl For low transmission, there is a huge 
increase of the current at-.small voltage leading to zero- 
bias conductance peaks uB 

Next, we study s/d xy junctions, see Figs. 2(b)-(c). 
As the s-wave gap A s increases, the s-wave gap more 
and more " shadows" the zero-energy MGS on the d-wave 



side. Since the width of the MGS decreases with de- 
creasing transmission T this effect is more pronounced 
for low transmission. Subharmonic gap structure (SGS) 
appears at voltages A s /n, n > integer, see Fig. 6(a) 
in Ref. [r]. All pthar SGS predicted to occur in asym- 
metric junctionsE_§E£l is partly washed out after angular 
averaging since it involves also the gap of the <i-wave su- 
perconductor (which depends on 6). The SGS at voltages 
A s /n depends on n being even or odd. For odd n the 
scattering states originating from the left (s-wave) su- 
perconductor will hit the resonant MGS at the <i-wave 
side, leading to negative differential resistance. For even 
n the scattering states responsible for SGS do not hit the 
MGS, and therefore this case resembles more what one 
finds in junctions with s-wave superconductors on both 
sides. When the transmission decreases, all SGS besides 
n = 1 vanishes, in which case the negative differential 
conductance is large. 

We also study the case when there are d-wave super- 
conductors on both sides of the junction. First, the case 
d x 2_ y 2 /d x 2_ y 2 is shown in Fig. |2|(d). In this case the 
MGS are absent since A (0) = A ((9). We can therefore 
compare our low-transmission curve [the case T = 0.026 
in Fig. 0(d)] with the tunneling calculation shown in 
Fig. 2 of Ref. [13], and we find that the curves are the same. 
At small voltage V there is a V 2 dependence of the cur- 
rent, and one finds a threshold voltage at eV = 2Ao. For 
intermediate transmission [the case T = 0.38 in Fig. @(d)] 
SGS appears at voltages 2Ao/n. The SGS is not as sharp 
as in the s-wave case because of angular averaging. 

For the d xy /d xy configuration, see Fig. ||(e), there are 
MGS's on either side of the junction. The difference com- 
pared to the d x 2_ y 2 / d x 2_ y 2 case is that the presence of 
MGS lower the threshold voltage to eV w A . This 
behavior is similar to the case when a BrcJ-UWigner res- 
onance is embedded in the normal region. ta The reason 
for the lower threshold voltage is that the MGS trans- 
mission resonances at the right (left) hand side are most 
pronounced for 9 w 7r/4. At this angle a quasipar- 
ticle coming from the left (right) experiences the gap 
A(7r/4) = Ao, resulting in a threshold voltage eV « Ao 
to make it possible for the quasiparticle from the left 
(right) side to hit the MGS on the right (left) side. Angu- 
lar averaging smears out the threshold. For small voltage 
V there is a V 2 dependence of the current in the low- 
transmission limit. For intermediate transmission SGS 
again shows up at 2A /n, somewhat washed out because 
of angular averaging. In the s/d xy case the negative dif- 
ferential conductance comes from resonant conduction of 
the right-moving scattering states through the MGS of 
the right-hand side. In the d xy /d xy case both left-movers 
and right-movers will pass through MGS (since there are 
MGS on both sides), and the resonant contribution from 
both left-movers and right-movers cancels. Therefore, 
there is in the d xy /d xy case no negative differential con- 
ductance. 

In Fig. ||(f) we show the I—V curves for the d x 2_ y 2Jd xy 
case. This case is similar to the s/d xy case shown in 



Figs. |(b)-(c) since Ai(0) = Ai(0). One important dif- 
ference is that the angular averaging now results in non- 
zero current for small voltages. Still there is negative 
differential conductance close to the voltage eV « Ao 
in the low-transmission limit. One could understand the 
I — V curve in Fig. |2|(f) as an averaging of the curves 
in Figs. 0(a)-(c). For small transmission the current is 
proportional to the voltage V in the small- voltage limit. 

In summary, we have calculated I — V curves for pla- 
nar (i-wave junctions. In some arrangements we find 
considerable negative differential conductance and lower 
threshold voltage for non-zero current due to resonant 
conduction through MGS. These features have not been 
discussed before in the literature and they should be ex- 
perimentally observable if superconducting gaps with d- 
wave symmetry exist in nature. 

I thank M. Samanta for many discussions on <i-wave 
superconductivity and for a lot of inspirational energy to 
carry this project to an end. Also, I have benefitted from 
discussions with P. F. Bagwell, Z. Ivanov, G. Johansson, 
V. S. Shumeiko, G. Wendin and S. Ostlund. 
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